SU-HEP-010 



S-waves and the Measurement of CP Violating Phases in B s Decays 

S. Stone and L. Zhang 

Physics Department, Syracuse University, Syracuse, New York 13244 

(Dated: December 18, 2008) 

Abstract 

Heavy, as yet undiscovered particles, can affect measurements of CP violation in the B system. 
Measuring CP violation in the B s system provides an excellent place to observe such effects since 
Standard Model sources are predicted to produce very small effects. The angle — 2/3 s , the "phase of 
B s — B s mixing," thought to be best measured in B s — > J /ip4> decays is of order -0.04, while the CP 
violating asymmetry in B s — > cfxf) is predicted to be zero, due to the cancellation of the mixing phase 
with the decay phase. Recent measurements of f3 s in J/ipcft, while not definitive, are much larger than 
the Standard Model predictions. Measurements in the B° and Df systems of analogous modes point 
toward a 5-10% contamination of S-wave K + K~ under the <f> peak. This S-wave was not taken into 
account in these recent analyses. Furthermore this S-wave can also materialize as a /o(980) meson 
that decays to tt + it~ , making the final state J/V'/o useful for measuring (3 S with the added advantage 
of not requiring an angular analysis. Rate estimates, while not precise, predict four to five times fewer 
such events than those in the J mode. The error on (3 S , however, may be similar. We also remark 
on S-wave problems with the B s — > <p<j) mode, and possible systematic checks using B s — > (frfo. 

PACS numbers: 13.20.He, 14.40.Nd 
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Measurements of Charge-Parity (CP) violation in the B meson system are sensitive to the 
presence of heavy, as yet undiscovered, particles. While CP violation has been studied exten- 
sively in the B° system where it is a large effect, very little information exists for the B s 
system. The CP violating angle —2(3 S , the so called "phase of B s — B s mixing" is a particularly 
important place to look for physics beyond the Standard Model, since the expected asymmetry 
is very small, sin(2/3 s ) = 0.037 ± 0.002 [3], thus allowing the effects of any new physics to be 
more easily observed. Both CDF ^ and DO 0] have investigated — 2(3 s using B s — > J/if)<f> 
decays. Central values have been found far from the expected Standard Model values, but the 
errors are large and the significance is in the 2-3<r range. 

Since the final state consists of two spin-1 particles, the final state is not a CP-eigenstate, 
yet it is well known that CP violation can be measured using angular analyses The analysis 
that have been heretofore carried out have ignored the possibility of an S-wave K + K~ system 
in the region of the 0. This S-wave can bias the result and not accounting for it certainly makes 
the error smaller. In fact, the analogous channel in B^ decay J/if)K*° is well known to have 
an S-wave Ktt component in the K* mass region. This interference, in fact, has been used by 
BaBar to measure cos(2/5) and thus remove an ambiguity in the value of (3 from the sin (2/3) 
measurement. The S-wave component in the region of the K* is measured as ~ 8% [6j. 

Perhaps it may be hoped that the S-wave K + K~ under the in J/if)<p ls smaller due to 
the relatively narrow width (r) of the (4.3 MeV) compared to the K* (51 MeV), but we will 
show here that this likely is not the case. (The S-wave problem was brought up at the recent 
CKM workshop in the discussion of the CDF and DO results Q.) This adds another amplitude 
and phase to the number of parameters that need to be determined in the full transversity 
analysis and increases the error over previous expectations [§[ . Similar considerations apply to 
the measurement of CP violation in the process B s — > 00. Here the problem is exacerbated by 
the presence of two 0's in the final state. The decay diagrams for both of these processes are 
shown in Fig. [TJ In both cases the ss forms a 0. Other manifestations of ss quarks are the rj, 
1]' and /o(980) mesons. The first two are pseudoscalars, while the last is a scalar. 




FIG. 1: Decay diagrams (a) B s — > J/ip<f>, and (b) B s — » 00. 



Final states with J /if) plus a spin-0 object have the J /if) fully polarized in the (J, J z ) (1,0) 
state, and thus an angular analysis is unnecessary. Unfortunately, the measurement of (3 S using 
the pseudoscalars has less sensitivity than using the at hadron colliders, because the large 
decay modes of the rj contain at least two photons and the rf at least one. Due to the relatively 
poor photon detection efficiency, larger backgrounds, and poor mass resolutions using photons 
at current experiments, the accuracy on (3 S is much poorer using the r\ or rf modes than in the 
mode, even though the angular analysis can be avoided. The scalar / state, however, hasn't 
been previously considered. 

We note that the /o(980) is an elusive object that decays largely into 7r + 7r~ but can also decay 
into K + K~. Some information about this particle can be gleaned from analysis of exclusive 
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.D+ final states. The 7r + 7r + 7r~ and K + K~7r + are of prime interest. The simple spectator decay 
diagram that results in the 0tt + or the /o 7r+ final states is shown in Fig. [2j We note that the 
ss system can form only isoscalar final states. 




FIG. 2: Decay diagram for a decay into a tt + plus an ss system. 

Of course, not only these final states are produced. A full Dalitz plot analysis, however, 
reveals the structure. The most recent and highest statistics study of the 7r + 7r + 7r~ final state 
has be done by BaBar [9]. The invariant mass spectrum and Dalitz plot are shown in Fig. [3j 
There best fit results showing fractions and phases are summarized in Table [H 



Decay Mode Decay fraction(%) 


Amplitude 


Phase(radians) 


/ 2 (1270)tt+ 
p(770)vr+ 
P (U50)tt + 
S-wave 


10.1 ±1.5 ±1.0 
1.8 ± 0.5 ± 1.0 
2.3 ± 0.8 ± 1.7 
83.0 ± 0.9 ± 1.9 


1. (fixed) 
0.19 ±0.02 ±0.12 
1.2 ±0.3 ±1.0 


0. (fixed) 
1.1 ±0.1 ±0.2 
4.1 ±0.2 ±0.5 


TOT. 

X 2 /NDF 


97.2 ±3.7 ±3.8 
437 _ 12 

422-64 — L '* 







TABLE I: Results from the BaBar — > ir + ir tt + Dalitz plot analysis. The table reports the fit 
fractions, amplitudes and phases. Errors are statistical and systematic, respectively. 

Thus, most of this final state is S-wave (83%). The mass projection of the Dalitz plot shown 
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FIG. 3: (a) The tt + tt + tt invariant mass distribution, (b) The symmetrized Df — > 7r + 7r + 7r Dalitz 
plot. From the BaBar collaboration 
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FIG. 4: Dalitz plot projections (dots with error bars) and fit results (solid histogram). (a) 
m 2 {'K + TT^)L ow , (b) rn 2 {TT + ir~)m g } l , (c) total m 2 {^ + ir~), (d) m 2 (7r + 7r + ). The hatched histograms 
show the background distribution. From the BaBar Collaboration [9(. 



in Fig. H] shows a relatively narrow peak in 7r + 7r mass within ~ ±90 MeV of the fo mass. This 
is in agreement with a previous analysis by FOCUS 



A Dalitz plot analyses of the K + K tt + final state have been performed by E687 [12j]. They 
find the results shown in Table [Til While extracting a precise ratio of decay rates is difficult, 



Decay Mode Decay fraction(%) Phase(degrees) 

K*(892)°K + 0.478 ± 0.046 ± 0.040 (fixed) 

<f)Tr + 0.396 ± 0.033 ± 0.047 178 ± 20 ± 24 

/o(980)7r+ 0.110 ±0.035 ±0.026 159 ± 22 ± 16 

/j(1710)vr + 0.034 ± 0.023 ± 0.035 110 ±20 ±17 

K*(U30)°K + 0.093 ± 0.032 ± 0.032 152 ± 40 ± 39 

Goodness of Fit -21n£=-1075 conf. level=80.2% 

TABLE II: Results from the E687 Df -> K+R-1T+ Dalitz plot analysis. The table reports the fit 
fractions, and phases. 
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FIG. 5: Dalitz plot projections for K + K invariant mass in — » K + K ir + from the CLEO 



Collaboration 13j. The signal is extracted individually in each mass bin, thus there is no background. 
The data are fit with a Breit-Wigner signal function for the 4> convoluted with a Gaussian for detector 
resolution and linear representation of an S-wave component (dashed line). The solid curve shows the 
sum. (Only the data is ascribed to CLEO, the fits have been added.) 



because the phases of the /o 7r+ and 07r + amplitudes are almost equal, within error, we can infer 
that the ratio is 

T(D+ -> f (980)7T + -> K + K-n + ) 

-±-^ ^ J - -^ = 0.28 ±0.12 . 1 

Clearly one cannot ignore the S-wave contribution here. However, this analysis is done over all 
of phase space and we concerned only with the low mass region. 

CLEO has looked explicitly at the low mass K + K~ region in the K + K~tt + Dalitz plot 
[13|. Fig. [5] shows the data in the region near 1 GeV. The signal is extracted by fitting the D s 
yield in each bin of K + K~, so no background remains in the plot. There is clearly an extra 
rather flat component of signal beneath the <fi. To estimate the size of this component, we 
have fit the CLEO data to a Breit-Wigner to describe the 0, convoluted with a Gaussian for 
detector resolution, and in addition added a linear S-wave component. The fraction of S-wave 
depends on the mass interval considered. For ±10 MeV around the mass we have 6.3% S-wave 
contribution, which rises to 8.9% for a ±15 MeV interval. (Note that these fractions depend 
on the experimental resolution.) 

We now try to compare the S-wave in D s decays with the one in B s —> J/ipK + K~ . The near 
equality of the fo and <j) masses removes this from consideration. Let us consider the energy. In 
the case of D s — > K + K~n + we have 1.97 GeV minus the pion mass, or 1.83 GeV available. In 
the case of B s we have 5.37 GeV minus 3.10 GeV or 2.27 GeV. A bit more energy but not that 
different. What does differ is that we have a 0~ pion in one case and a T J/tp in the other 
case, and the spin may matter. 

Now we wish to predict how many events we may get should we look for B s — ► J/tpfo with 
fo -> 7T+7T-. CLEO-c measures B(D+ -» K + K~-k + ) = (5.50 ± 0.23 ± 0.16)% 0- First we 
need to derive a <f)7r + branching ratio (with — > K + K~). CLEO-c gives a range of branching 
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ratios computed for varying cuts around the mass. Extrapolating their data to a zero width 
interval allows us to remove the S-wave component under the (see Fig. [5]). This gives a value of 
B(D+ — > 07r + , —> K + K~)—(l.d±0.1)%. We also use the CLEO absolute branching fraction 
measurement B(Df — > 7r + 7r~7r + ) = (1.11 ± 0.07 ± 0.04)%, and the BaBar Daltiz analysis of 
the 37r mode to derive a branching fraction for B(D+ — > /o 7r+ ~~ ► vr + 7r + 7r _ ) = 0.30%, where 
we found by integrating that (27±2)% of the 3n mode is due to a narrow f within ±0.13 
GeV 2 of the fo mass-squared (~ ±90 MeV of the fo mass of 980 MeV). This comes from 
an examination of the BaBar data shown in Fig. SJc). Thus, we expect a sample of events 
B s — > J/ipf , /o — > 7r + 7r~ that is (19±2)% of the sample. If we include the contribution 

of all the 7r + 7r~ S-wave combinations above background in the three peak bins, our estimate 
would increase to (27±3)%, as there is significant 7r + 7r~ S-wave outside of the fo mass region. 

We now make another estimate of the relative fo — > 7r + 7r~/0 — > K + K~ rate. Semileptonic 
.D s decays provide another basis of estimate. Here the ss pair is produced opposite a virtual 
W~ , and thus is produced opposite an object with the same spin as the J/if). Measuring the 
ratio of semileptonic decay rates where the final state meson has the maximum momentum (at 
q 2 = 0, where q 2 is the 4-momentum transfer) would approximate the available energy quite 
well (see above). CLEO has recently measured the branching ratios integrated overall q 2 [14| 

B(Dt -> 0e + i/), -> K + K- = (1.1 ± 0.2)% (2) 
B(D+ -> / e + z/), /o-Tr+TT = (0.13 ± 0.04)% . 

(3) 

Unfortunately, the q 2 distributions are not yet available so we will use the integrated rates; we 
know the maximum rate is at q 2 of zero. The ratio of observable fo/4> here is (13±4)%, where 
the error is only from the CLEO data and not the model approximations. This estimate is very 
qualitative in that we do not have the q 2 equals zero data, but it is consistent with our previous 
estimate, given the limitations in the data. 

Thus, using our value from hadronic D s decays, we expect that 

r(fi s -> -> K+ir-) 1 ' 

This rate may increase to 30% if there is additional S-wave under the fo. We note that 
any 7r + 7r~ S-wave, or even D-wave, does not constitute a background for the CP violation 
measurement as it has the same CP as the fo; also, since the o° doesn't have any ss in its 
wave-function, it is unlikely that the ss system forms a 1~ state la ]. The purity of the sample 
can be checked by examining the polarization of the leptons from the J/ip — > £ + £~ decay; this 
serves as an experimental check on the amount of opposite CP P-wave. For pure fo, in the limit 
of massless leptons the angle 9 of the l + in the J / ip rest frame with respect to the fo direction 
in this frame, must be distributed as sin 2 9. 

In the decay B s — > <fxp, the decay phase cancels the mixing phase if only Standard Model 
particles are involved, so this reaction provides an interesting place to look for manifestations 
of new physics [16]. The presence of the S-wave in the K + K~ distributions at low mass also 
applies to the 00 final state. Here, however, there are two possible S-wave contributions that 
must be taken into account since there are two 0's. It may be that using the 0/o final state 
can help in understanding this effect. Here an angular analysis is only necessary to disentangle 
the amount of S-wave under the peak. We note however, that events with two charged kaons 
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and two charged pions can also result from the K K*° final state, so a great deal of care must 
be exercised. (These final states seem to be well separated by the large phase space available, 
at least at first glance.) 

In conclusion, we predict on the order of a 10% K + K~ S-wave contribution that contami- 
nates B® — > J/ip(f), <t> — > K + K~ and must be taken into account by adding an additional S-wave 
amplitude and phase to the parameterization of the decay width. Inclusion of this S-wave will 
change the central values and increase the errors in current analyses of f3 s . We suggest a new 
decay B s — > J/ipfo, fo — > tt + it~. Rough estimates indicate that we can expect at least ~l/5 
as many events as in J/-00. An angular analysis is not necessary to find f3 s using these events, 
and thus the statistical error on j3 s could be comparable to the and the systematic error smaller 
than those for the J/ip<f) mode. A critical look needs to be given towards the angular analysis 
of the B s — ■> 00 mode since the S-wave can enter twice. It may be possible to access the same 
physics using the B s — > 0/ o , where an angular analysis may only be necessary to ascertain the 
fraction of S-wave under the peak. Measurement of CP violation in the B s system continues 
to be of utmost importance. 

We would like to thank M. Artuso, J. Rosner and several LHCb colleagues for useful discus- 
sions. We thanks the U. S. National Science Foundation for support. 
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